Kondo-lattice theory of anisotropic singlet superconductivity 
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Kondo-lattice theory for the t-J model and a phenomenological theory based on it are developed to 
study superconductivity in the vicinity of the Mott metal-insulator transition. Since the quenching 
of magnetic moments by single-site quantum spin fluctuations or the Kondo effect is reduced by the 
opening of a superconducting gap, spin density wave (SDW) can appear in a superconducting state 
and the Knight shift can deviate from the Yosida function to be small. The electron-phonon interac- 
tion, which arises from the modulation of the superexchange interaction by phonons, is crucial in the 
coexistence of superconductivity and SDW. It is proposed that the coexistence of superconductivity, 
a double-Q SDW, a double-Q lattice distortion, and a double-Q charge density wave induced by 
the SDW rather than the lattice distortion, is responsible for the checkerboard structure and the 
zero-temperature pseudo-gap observed in under-doped cuprate superconductors. 

PACS numbers: 74.20.-z, 71.10.-w, 71.30.+h 



I. INTRODUCTION 

High temperature (high-T c ) superconductivity is one of 
the most interesting and important issues in condensed- 
matter physicSfi*2i2ii5i£ It occurs in highly anisotropic 
quasi-two dimensional cuprate oxides, i.e., on Cu02 
planes. Parent cuprates are Mott insulators. Antifcr- 
romagnetism occurs at low temperatures because of the 
superexchange interaction between nearest neighbor Cu 
ions on Cu02 planes, which is as strong as 7 - 



J= -(0.10-0.15) eV. 



(1.1) 



When holes or electrons are doped, cuprates become met- 



als and eL 



-wave or <i7-wave superconductivity occurs 



on CuC>2 planes. High-T c superconductivity occurs in the 
vicinity of the Mott metal-insulator (M-I) transition. 

The Hubbard model is one of the simplest effective 
Hamiltonians for the Mott M-I transition. On the basis 
of Hubbard's^ and Gutzwiller' o 10 ' 11 ! 12 theories, which 
are under the single-site approximation (SSA), it is spec- 
ulated that the density of state (DOS) has a three-peak 
structure, with the Gutzwiller band between the up- 
per Hubbard band (UHB) and the lower Hubbard band 
(LHB). Another SSA theory confirms this speculation, 13 
showing that the Gutzwiller band appears at the top 
of LHB for n < 1, with n being the electron density 
per unit cell. The supreme SSA, which considers all the 
single-site terms, is reduced to determining and solving 
self-consistently the Anderson model ) 14 ! 15 ' 16 which is an 
effective Hamiltonian for the Kondo effect. The three- 
peak structure corresponds to the Kondo peak between 
two subpeaks in the Anderson model. The Kondo effect 
has relevance to the Mott M-I transition. 

The supreme SSA is rigorous for infinite dimensions 
within the restricted Hilbert subspace where no order 
parameter exists^ It is later shown that the supreme 
SSA can also be formulated as the dynamical mean-field 
theoryi&i&22«2I (DMFT) and the dynamical coherent po- 
tential approximation^ 2 - (DCPA); they are exactly the 
same approximation as each other. 



As a function of doping concentrations S = 1 1 — n| , 
superconducting (SC) critical temperatures T c show a 
maximum at an optimal one 6 ~ 0.15, which is denoted 
by 5 op . Cuprates with S < S op , S ~ 5 op , and 8 > <5 op are 
called under-doped, optimal-doped, and over-doped ones, 
respectively. Specific heat coefficients of optimal-doped 
ones are as large a a 23 i 24 ' 25 



7 = 10-15 mJ/K 2 • Cu0 2 -mol, 



(1.2) 



which are consistent with the prediction of Gutzwiller's 
theor y 10 ' 11 ! 12 From Eq. (|1.2[) . the bandwidth of quasi- 
particles is estimated to be 



W* = 0.3-0.5 eV. 



(1.3) 



Cooper pairs can also be bound by a magnetic exchange 
interaction^ 6 - One of the simplest effective Hamiltonians 
for cuprates is the t-J model; 27 it is derived from the 
Hubbard model^ 6 - and is also derived from the the d-p 
model~ which considers 3d orbits on Cu ions and 2p 
orbits on O ions. It is quite easy to reproduce observed 
T c if the phenomenological J and W* are used. It is 
proposed by an early theor y 29 ! 30 that the condensation 
of G?7-wave Cooper pairs between the Gutzwiller quasi- 
particles is responsible for high-T c superconductivity. 

There are pieces of evidence that the electron-phonon 
interaction is strong on Cu02 planes: the softening of the 
half-breathing modes around (±7r/a, 0) and (0, ±7r/a) in 
the two-dimensional Brillouin zone ^D-BZ), 3 ^! 3 ^ 3 .! 3 ^ 5 . 
with a the lattice constant of the Cu02 planes, the soft- 
ening of Cu-0 bond stretching modes around (±7r/2a, 0) 
and (0, ±7r/2a) in 2D-BZ ; 36 ' 37 kinks in the dispersion re- 
lation of quasi-particleS ) 38 ' 39 and so on. On the other 
hand, observed isotope shifts of T c are small^ - which 
implies that the strong electron-phonon interaction can 
play no significant role in high-T c superconductivity it- 
self. 

The modulation of the superexchange interaction by 
phonons gives a strong electron-phonon interaction ! 41 ! 42 
It can explain the softening of phonons. The virtual ex- 
change of a relevant phonon gives a mutual attractive 
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interaction. Since it is strong between next nearest neigh- 
bors but is weak between nearest neighbors, it can play 
no significant role in the binding of cfy-wave Cooper pairs, 
which arc bound mainly between nearest neighbors. 

Many anomalies are observed in under-doped 
cuprates; some of the most crucial ones are the 
pseudo-gap,— ' 43 ' 44 ' 45 ' 46 ' 47 ' 48 the stripe or checkerboard 
struct ure; l i 6 i 49 i 5[l i 51 i 52 i 53 i 54 the zero-temperature pseudo- 
gap (ZTPG); 6 -! 5 ^ 5 ^ 6 - and so on. One may argue that 
an exotic mechanism must be responsible for not only 
the anomalies but also superconductivity. Since T c 's are 
more or less low when one or some of the anomalies are 
observed, the other may argue that the role of what 
cause the anomalies is simply the suppression of T c . The 
anomalies should be consistently explained with high-T c 
superconductivity. 

The supreme SSA is mapped to solving the Anderson 
model, as is discussed above. When no SC gap opens in 
the Hubbard, t-J, or d-p model, the Fermi surface (FS) 
exists in the Anderson model so that the ground state of 
the Anderson model is the normal Fermi liquid (FL). In 
the SSA, the ground state is stabilized by the Kondo ef- 
fect and is the FL. A perturbative treatment of intersite 
effects starting from the non-perturbative SSA theory is 
simply Kondo-lattice theory. Beyond the SSA, the FL is 
further stabilized by the Fock-type term of the superex- 
change interaction ) 57 ' 58 which favors a local singlet or a 
resonating valence bond (RVB) on each pair of nearest 
neighbors; the stabilization mechanism is physically the 
same as the RVB mechanism . 27 ' 59 ' 60 

Critical fluctuations make T c down to +0 K in two 
dimensions, 61 which implies that they are also cru- 
cial in highly anisotropic quasi-two dimensions. If the 
anisotropy is large enough but T c is still high enough, SC 
critical fluctuations can cause the opening of a pseudo- 
gap in the SC critical region i 62 ' 63 

The Kondo effect is the quenching effect of magnetic 
moments by single-site quantum spin fluctuations, whose 
energy scale is the Kondo temperature fcsTk or Tk. 
When an SC gap opens, the FS vanishes in the Anderson 
model so that Tk is reduced. It is possible that Tk is 
so reduced in an SC state that spin density wave (SDW) 
can appear in the SC state. One of the purposes of this 
paper is to develop microscopic and phcnomenological 
theories to study possible coexistence of superconductiv- 
ity and SDW in under-doped cuprate superconductors. 
The organization of this paper is as follows: In Sec. HH 
Kondo-lattice theory of superconductivity is extended 
such that it can treat properly the reduction of Tk in 
an anisotropic SC state. It is shown in Sec. IIIII that, be- 
cause of the reduction of Tk , SDW can appear in the SC 
state and the Knight shift can deviate from the Yosida 
function. In Sec. IIV( a phenomenological theory is de- 
veloped on the basis of the microscopic theory developed 
in Sec. [TT] to study the coexistence of superconductivity 
and SDW; it is shown that the electron-phonon interac- 
tion is crucial in the coexistence. It is proposed in Sec.lVl 
that the coexistence is responsible for the checkerboard 
and ZTPG observed in under-doped cuprate supercon- 



ductors. Discussion is given in Sec. IVIl Conclusion is 
given in Sec. IVIII Anderson's compensation theorem is 
studied in Appendix lAl 

II. KONDO-LATTICE THEORY 

A. Nambu representation 

In this paper, the t-J or t-J-Uoo model on the square 
lattice is considered: 

icr i=jLj a (ij) 

i 

with jiia = a ia ai ai e a the band center, and 

Si = Y,\i°f^f>°f) a l^ (2-2) 

a/3 

with a x , tjy, and a z being the Pauli matrixes. The dis- 
persion relation of electrons is given by 

Ma 

with N the number of unit cells. The bandwidth of E(k) 
is W = (3-4) eV for cuprates. The third term in TL is 
the superexchange interaction between nearest neighbors 
(ij). According to field theory, it arises from the virtual 
exchange of a pair excitation of electrons between UHB 
and LHB^ it is certain that the superexchange inter- 
action is present in not only insulating but also metallic 
phases. In the Hubbard model, the well-known value of 
J = —4:t 2 /U, with t the transfer integral between nearest 
neighbors and U the on-site repulsion, can only be repro- 
duced if the bandwidths of UHB and LHB are ignored. 
When actual non-zero bandwidths are considered, \J\ is 
much smaller than 4t 2 /U . The experimental J can only 
be explained beyond the Hubbard model, for example, 
by the d-p model; 6 - 5 - In this paper, Eq. ([1.1)1 is assumed 
for J. It is convenient to define 

J s (q) = -2 J [cos(q x a) + cos(q y a)} , (2.4) 

with J = —(0.10-0.15) eV and a the lattice constant. 
In the last term in 7i, Uoc is an infinitely large on-site 
repulsion such as Uoo/W — > +oo, which is introduced 
to exclude any double occupancy. Weak inter-layer cou- 
plings or quasi-two dimensionality is implicitly assumed 
so that SC T c can be non-zero. 

It is definite that, when n ~ 1, T c of dj wave is much 
higher than T c of other waves i 29 ' 30 Provided that only a 
d7-wave SC order parameter exists, which is denoted by 
Ad 7 (k), the single-particle Green function in the Nambu 
representation is rigorously given by 
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&(ien,k) = - J^' kBT dre^ (t t (j^^j (4»> a-k-«r(r))\ 



ie n - S(k) + /i - S CT (i£„, k) 
-AS 7 (k) 



-A rf7 (k) 



ie„ + J5(k) - fi + Y,- a {-ie n , k) 



(2.5) 



with a^pij) and ako-(T) being creation and annihilation 
operators in the Heisenberg representation, /i the chem- 
ical potential, and E CT (ie„,k) the self-energy. The deter- 
minant of Q~ 1 (ie n , k) is given by 

\g~ x (ie n ,k)\ = [i£„-.E(k)+/i-£ (7 (i£ n ,k)] 

x [ie n + E(k) -fi + S_ CT (-i£ n , k)] 



It follows that 

Ga(iSn, k) = 

with 

G CT (ie„,k) = 



|Ad 7 (k)| 



G CT (ie„,k) F(«£ n ,k) 
F*(ze n ,k) -G_ <r (-ie n ,k) 



ze n + g(k) - (j, + S- tT (-ie K , k) 



and 



F(ie n ,\s) 



A rf7 (k) 
1 (ie n , k)| 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



B. Mapping to the Anderson model 

The Feynman diagrams of the t-J model are classified 
into single-site and multi-site ones. If only Uoo and the 
single-site Green function, which is defined by 



Ra{ie n ) = —y)G a (ie n ,k), 



(2.10) 



appear in a diagram, it is a single-site one. If an inter-site 
part of G,j(«£ n ,k), F(«£„,k), or J appears in a diagram, 
it is a multi-site one. According to this classification, the 
self-energy is divided into single-site and multi-site ones: 



T, a (ie n , k) = E a (is n ) + AT, a (ie n , k). 



(2.11) 



The single-site self-energy 'E (7 (ie n ) is given by that of the 
Anderson model: 



a = td E nd ° 



kcr 



Fc(k)c ka .c kff + Uoon d ^n dl 



1 



E M 



k^ + V k d W* 



(2.12) 



with rido = d^da, td the energy level of d electrons, E c (k.) 
the dispersion relation of conduction electrons, N the 
number of lattice sites, Vk the hybridization matrix be- 
tween d and conduction electrons, and the on-site Uqo 
being the same as that for the t-J model. The Anderson 
model is characterized by e d — fx, with /2 the chemical 
potential, and the hybridization energy defined by 

A ( £ ) = iEl^| 2 <^ + M-£c(k)]. (2.13) 

k 

When A(0) > 0, the FS exists in the conduction band 
E c (k). Then, the ground state of the Anderson model 
is the FL at least when A(0) > 0. Consider the Green 
function defined by 



1 



M - e d ~ S ff (i£„) - 3 fde' . ^ £ — 



16 n 



(2.14) 

Provided that — fl and A(£) are determined to satisfy 



RaiiSn) — G CT (i£ n ), 



(2.15) 



the single-site T, a (ie n ) is simply the self-energy of the An- 
derson model. Equation (|2.15p is the mapping condition 
to the Anderson model; the multi-site AE CT (i£ ra ,k) and 
Ad 7 (k) should also be self-consistently calculated with 
the single-site T, a (ie n ) to satisfy Eq. (|2.15p . 

When AE CT (j£ n , k) and A^ 7 (k) are ignored in the map- 
ping condition (|2.15[) . the approximation is simply the 
supreme SSA, the DMFT, or the DCPA, which considers 
all the single-site terms. In this paper, both AS cr (?£ n ,k) 
and A(2 7 (k) are considered to treat the SC state properly. 

According to Eq. (|2. 15|) . DOS of the t-J model is the 
same as that of the Anderson model: 



P{e) 



- Irni^e + iO) 



-ImG CT (£ + iO). (2.16) 



The electron density of the t-J model is also the same 
as that of the Anderson model: n — (l/N) ^2 irJ (n ilT ) = 
y] (rido-)- It follows from the mapping condition (|2.15|) 
that 



£d - M = ( a - Ml 



(2.17) 



and 



^ S . (2-18) 
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When T = K, Im£ CT (-HO) = so that 



A(0) 



7Tp(0) 



[ReR a {+iO)} 2 + [np(0)Y 



(2.19) 



When no SC gap opens, p(0) > so that A(0) > and 
the ground state of the Anderson model is the FL. When 
an SC gap opens, p(0) = so that A(0) = or no FS 
exists. When the SC gap is anisotropic, the ground state 
of the Anderson model is still the FL, as is examined 
in Sec. IIII Al The Kondo temperature 2k or fcs2k is a 
characteristic low-energy scale, which is precisely defined 
in Sec. Ill CI It is assumed in the following part of this 
paper that T <C 2k and T c <C 2k. In the presence of an 
infinitcsimally small Zceman energy —ah in the Anderson 
model, the self-energy of the FL is expanded as 

+0(e 2 /k B T K ) + 0(k B T 2 /T K ) , (2.20) 

with 7 and (f> g all being real. When n ~ 1, charge 

fluctuations are suppressed so that <^) 7 3> 1, S S> 1, and 



<t>s/<i>-> 



(2.21) 



i.e., the so called Wilson ratio is about 2.— 

When AS IT (ie n ,k) is considered but A £ j 7 (k) is ignored 
in Eq. (|2.7p , the Green function is given by 



G a {ie n ,\z) 



1 



1 



i 7 ie n - £ (k) - AS CT (ie„, k) 
+ [incoherent term], (2.22) 



with 



and 



£o(k) = j-[E + £(k)-/i], 



AE CT (i£„,k) = J-A£ CT (ie„,k) 



(2.23) 



(2.24) 



The first term in Eq. (|2.22p is the coherent term, 
which describes the Gutzwiller band with band-width 
W* = 0(k B T K ) at the top of LHB. The second term in 
Eq. (|2.22p is the incoherent term, which describes LHB; 
UHB is infinitely high in the t-J-Uoo model. 
It is convenient to define 



A d7 (k) = J-A d7 (k) 



(2.25) 



and 



D a (ie n ,k) = [ie n - £ (k) - A£ .(ie n ,k)] 

x [fe„ + £ (k) + A£_ CT (-ie n ,k)] 
-|A rf7 (k)| 2 . (2.26) 



When both of AS -(ze„,k) and A^ 7 (k) are considered in 
Eq. (|2.7[) . the Green function is given by 



with 



and 



G CT (ze„, k) = , (2.28) 

A d7 (k) 



F(ie n ,k) 



D a {is n ,\s) ' 



(2.29) 



The incoherent part is ignored here and in the following 
part. 



C. Intersite exchange interactions 

The polarization function in spin channels is also di- 
vided into the single-site and multi-site ones: 



7r s (iuji , q) = 7T S (iuji ) + Att s (iui , q) . 



(2.30) 



The spin susceptibility of the Anderson model is given 
by 



2TT s (iu)l) 



1 - UooTTs^itOl) ' 

and that of the t-J model is given by 

27r s (iw;,q) 



Xs(*Wi,q) 



1 - [\J.M) + ^oo] 7T s («w ; , q; 



(2.31) 



■, (2.32) 



with J s (q) defined by Eq. (|2.4p . It should be noted that 
n^iuJt) = 1/E/oo + O [l/Ulxs(i"i)] , (2.33) 

and 

Air s (iu h q) = O [1/Ulxs(m, q)] . (2.34) 

A physical picture for Kondo lattices is that local spin 
fluctuations on different sites interact with each other 
by an intersite exchange interaction. Then, an intersite 
exchange interaction I s {iiOi, q) is defined by 



X*(«wj,q) 



1 - jl s (iuj[,q)x s (iu)i) 



(2.35) 



It fol lows from Eqs. (j2~3T|) . (|2~321) . (jOg)) . (|2~33|> . and 
(|235|) that 

/.(iwj.q) = J s (q) + 2C/^A7r s (z Wi ,q), (2.36) 
in the limit of Uqo/W — > +oo. 
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The Kondo temperature is denned by 
fc B T K = r lirn K [l/x s (0)]. 



(2.37) 



Since the mapped Anderson model is determined for a 
given T, it includes the given T as a parameter in addition 
to the temperature T itself. In Eq. (|2.37[) . the limit of 
T — > K stands for the limit for the temperature T; 
the parameter T should be the given T. The Kondo 
temperature Tk depends on the parameter T. 

The main term of 2U% <D Air s (iuji, q) in Eq. (|2.36p is an 
exchange interaction arising from the virtual exchange 
of a pair excitation of Bogoliubov's quasi-particles in the 
SC state. When the reducible and irreducible three-point 
single-site vertex functions in spin channels are denoted 
by A s (is n + iuJi,ie n ;iLUi) and X s (ie n + ioj u ie n ;ibJi), re- 
spectively, it follows that 



A s (ie r , 



\ s (ie n +iuji,i£ n ;iuji) 

1 - UooTT s (ilUl) 



(2.38) 



according to the Ward relation^ Since S — A s (0, 0; 0), 
it follows that 



X s {i£ n + iu>i,i£ n ;iwi) 



20s 



(2.39) 



for £ n — y and uji — » 0. When Eq. (|2.39|) is approxi- 
mately used, the exchange interaction is given by 



Jct(m,q) 



P(iui,q)-P(iuJi) , (2.40) 



with 



P(iu>i,q) = -k B T Jj [ G °-( i£ » + %u>u k + q)G ff (i£„, k) + F(fe„ + iwj, k + q)F(i£ n , k)] , (2.41) 

I 



which is derived in the random-phase approxima- 
tion (RPA) for pair excitations of Bogoliubov's quasi- 
particles. In Eq. (|2.40[) , the single-site term, 

(i£ n + iu>i)R a (i£ n ), (2.42) 

with 

Mis n ) = jjYl G*{ie n , k), (2.43) 

k 

is subtracted. When the so called mode-mode coupling 
term among various types of fluctuations, which is de- 
noted by — 4A(icL>;, q), is included, it follows that 

I s (iui,q) = J s (q) + J Q (iui,q)-4A(iwi,q). (2.44) 

The mode-mode coupling term corresponds to that in the 
so called self-consistent renormalization (SCR) theory of 
spin fluctuations ) 68 ' 69 which is relevant in the weakly cou- 
pling regime. 

When Eq. (|2.39[) is approximately used, the mutual 
interaction mediated by spin fluctuations is given by 

^(UX s ) 2 [ X s(t^,q) -Xs(m)} = \fii:(iu>i,q), (2.45) 
with A s standing for A s (0, 0; 0) and 

r*f \ I s (iui,q) , , 

/ s (iui, q) = r— ■ (2.46) 

1 - - l I s [iUJuq)Xs^i) 



In Eq. (|2.45|) , the single-site term is subtracted and two 
<j) s appear as effective three-point vertex functions. In the 
vicinity of the Mott M-I transition, the spin-fluctuation- 
mcdiatcd interaction is simply the exchange interaction 

D. Multi-site self-energy 

The multi-site self-energy AS (T (e + iO, k) is perturba- 
tively calculated in terms of I s (iu>i, q) or I*(iu>i, q). First 
of all, the Fock-type term of the superexchange interac- 
tion J s (q) should be considered to stabilize the normal 
or unperturbed statei 57 ' 58 When J s (q) is treated in the 
mean-field approximation, three types of order parame- 
ters are possible such as superconductivity, antiferromag- 
netism, and bond order (BO). This implies that antifer- 
romagnetism and BO compete with superconductivity. 
Not only SC fluctuations but also antiferromagnetic and 
BO fluctuations should be considered in calculating the 
multi-site self-energy AS CT (e + iO, k). 

When T = K or T is so low that intersite thermal 
fluctuations can never be developed, AS CT (e + iO, k) is 
expanded such that 

AY, a {£ + i0, k) = J- [AS (k) - A0 7 (k)e + •••]. (2.47) 

The dispersion relation of quasi-particles is given by 

£(k)=[£*(k)-/x]/^(k), (2.48) 
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with 



and 



E* (k) = S + E{k) + AE (k) , (2.49) 



,(k) = <£ 7 + A0 7 (k). 



(2.50) 



According to the FL relation, the specific heat at T c < 
T <C Tk is approximately given by 



C(T) = 7 T - 



(2.51) 



with 



f^4^£<5[M-£*(k)]0 7 (k) 



= =;n 2 k 2 B lP (O). 



Here. 



£ k %-.E*(k)]^(k) 
£ k %-£*(k)] 



(2.52) 



(2.53) 



is an average of 7 (k) over the FS and 

p(0) = ^£%-£*(k)]. (2.54) 



When the development of intersite quantum fluctuations 
is as large as that of local quantum spin fluctuations, it 
is possible that 



2<^ 7 



(2.55) 



It is also possible that this relation holds at not only 
T > T c but also T < T c . 

When intersite thermal fluctuations are much devel- 
oped, the expansion (|2.47[) is never relevant; precise T 
and e dependences of AT, a (ie n ,k) should be considered. 
Normal-state anomalies at T > T c such as those ob- 
served in under-doped cuprate superconductors can only 
be explained in terms of deviations from the expansion 
(|2.47p . For example, a large pseudo-gap opens provided 
that — ImAI] cr (e + i0,k) has a large and sharp peak at 
e = 0. 6 ^ 



E. Gap equation 



An average of /* {uj + iO, q) given by Eq. (|2.46p over a 
low-energy region \lo\ < /csTk is expanded as 



Bz(l*(u + i0,q)) ui =J2l(i,m)e 



i{q x l-\-q y 'm)a 



(2.56) 



i ni 



with (•••)„ standing for the average and I and m being 
integers. When the nearest-neighbor component, 



'(±1,0) — 1 {Q,±l) 



J (o,±i) = J l I 



(2.57) 



is only considered as the attractive interaction, two 
types of singlet superconductivity is possible such as 
anisotropic s and e?7 waves. In this paper, c?7-wave su- 
perconductivity is considered since T c of drf wave is much 
higher than T c of s wave. When ec(T) is defined by 



with 



A d7 (k) = - %7 (q) eG (T), (2.58) 



7y d7 (q) = cos^o) - cos(k y a), (2.59) 



being the form factor of cfy-wave Cooper pairs, the gap 
equation is simply given by 



i = ^ri(^) 2 ^E 



D a (ie n ,k) ' 



(2.60) 



with D a (ie n ,k) defined by Eq. (|2~26|l . where Eq. (J27 



should be used as A d7 (k). When the gap equation is 
self-consistently solved, the Green function is eventually 
given by Eq. (|2~27)) . 

In the Kondo-lattice theory developed in this paper, 
in principle, single-site properties such as £o(k), ^> 7 , and 
(f> s , the multi-site self-energy AE cr (ze„,k), the intersite 
exchange interaction Jj" between nearest neighbors, and 
the gap function ec (T) should be self-consistently calcu- 
lated with each other. The framework is in parallel to 
that of the conventional theory of superconductivity ex- 
cept that the unperturbed state or the single-site proper- 
ties should be self-consistently calculated in the supreme 
SSA with the multi-site or intersite properties. 



III. REDUCTION OF THE KONDO 
TEMPERATURE 

A. Spontaneous magnetizations 



In general, ReR a {+iO) ^ 0. According to Eq. (|2~T7)) . 
when the SC gap opens below T c , A(e) has also a gap 
structure; A(0) = at T = K. Although A(e) should 
be self-consistently determined to satisfy the mapping 
condition (|2. 15|) or (|2.17p . the gap structure of A(e) is 
phenomenologically treated. Since A(e) oc |e| for small 
|e| at T = K for c?7-wave superconductivity, 

Ax, -D<s<-e 

Ao + (Ai-Ao)|e|/|eo|, -e < e < 

A(e) = ^ Ao + (A 2 -Ao)|e|/|eo|, < e < e , 

A 2 , e < £ < D 

0, \e\ > D 

(3.1) 

is assumed; eo — (l/2)ec(0) and D ~ W. According 
to Eq. ([2~T7|) . Ai ~ 0{W/n) and Ai > A 2 , since the 
chemical potential is at the top of LHB. When the SC 



gap opens, A = at T = K. However, A is treated 
as another parameter. It is assumed that Ax > A > 0. 

Since the on-site Uoo is infinitely large, no doubly oc- 
cupied configuration of d electrons appears in any eigen- 
state. Then, the lowest singlet state of the Anderson 
model is expanded such thai^S 



®s = [A + ^2A d .. ka dic ka 



^ ^ ^kCT;po-'C ktr C po ./ 
ker per' 



|0>, (3.2) 



with |0) being the Fermi vacuum for conduction electrons 
with no d electron. It satisfies 

(n - ftrt + Wext)$., = e s $ s , (3.3) 

with H defined by Eq. (|27l2j) . 

£4.^, (3.4) 



N = 2 n do 




ktr 



FIG. 1: l/7rAiXs(0) as a function of n; (i) A /Ai = 0, (ii) 
Ao/Ai = 0.1, and (iii) A /Ai = 0.2. Solid, dotted, dashed, 
chain, and chain-dotted lines are for eo/Ai = 0, 0.1, 0.2, 0.3, 
and 0.4, respectively; every solid line shows Eq. ()3.13J| . In any 
line, l/7rAiXs(0) > even for n ~ 1. 



(3.5) 



being infinitesimally small external fields, and E s the en- 
ergy of the singlet. When only A^ka, and A kavpa f 
are considered, it follows thafcZi 



2 r 

E s = / deA(e) 



Ed - A^ - e — Eg 
2 ~' ( 3 - 6 ) 



(E d -Aii-e- E„y - h 2 
where Ed is the energy of the lowest doublet; 

A 2 , D 

E d ~ E + e d - \i In — , 

7T Al 



(3.7) 



with Eq the energy of the Fermi vacuum. 

When A > 0, it is obvious that E s < Ed- Even 
if Ao = 0, E s < Ed- The ground state of the mapped 
Anderson model is a singlet within the phenomenological 
model (|3.1[) . i.e., when A(e) oc |e| for small |e|. 

The d electron density is given by 

9E ° (3.8) 



It follows that 
n 



dAfi 



2 f« ^ A(e) 



1 - n nJ_ D (E d - e - E s y ' 



(3.9) 



which gives E s as a function of n. The magnetization of 
d electrons is given by 



dE s 

so that the susceptibility is given by 

^ (0) = hi 

4(1 -n) f° ^ A(e) 



(3.10) 



d [E d ~e- E a f 



(3.11) 



First, consider the case where no gap structure is de- 
veloped in A(e), i.e., Ao = Ai or eo = 0. It follows 
that 

E s =E d -^—<E d , (3.12) 



and 



X.(0) = 



2Ai 1 - n 



(3.13) 



Here, D — > +oo is assumed since D ^> \E d — E s \. 

When the gap structure is developed in A(e) or when 
Ao/Ai is small and eq/Ai is large, Xs(0) is enhanced 
from Eq. (|3.13|) or l/x s (0) is reduced. Figure [T] shows 
1/7tAiXs(0) as a function of n for various Ao/Ai and 
£o/Ai, where D — > +oo is also assumed. 

The Neel temperature Tn is defined by 



1 



1 



1 



J s (q)--J Q (0,q) + A(0,q) 



= 0, 



_X,(0) 4^ Jt=Tn 

(3.14) 

where Tn should be maximized as a function of q; the 
q that gives the maximal Tn is the wave number Q of 
antifcrromagnetic moments. If Eq. (|3.14[) gives a positive 
Tn, antiferromagnetism coexists with superconductivity. 

According to Fig. []] l/x s (0) is largely reduced when 
a large gap opens, for example, when Aq/A\ ~ 0, 
eo/Ai > 0.2, and n > 0.9. When l/x s (0) is reduced, 
Jq(0, q) and A(0, q) are also reduced. On the other hand, 
the superexchange interaction J s (q) is scarcely reduced. 
Provided that the gap is as large as cg(0)/Ai > 0.4, 
therefore, it is probable that antiferromagnetism appears 
for n > 0.9, at least, at T — K or antiferromagnetism 
coexists with superconductivity for n > 0.9. When J s (q) 
is given by Eq. (|2.4p , the wave number of antiferromag- 
netic moments is Q ~ (±7r/a, ±7r/o). 
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B. Small Knight shift 

The Knight shift is proportional to the static homo- 
geneous susceptibility X s (0,qo)j with q standing for an 
infinitesimally small q such as |q| — > 0. It follows that 



Xs (0,q o ;T) = 



1 



X (q )+X 1 (q ;T)' 



(3.15) 



with 

Y (q ) = 
and 

Xi(q ;T) 



J s (qo) 



/07)2 F(O) + A(O,q o ), (3.16) 



x2(o) 



(3.17) 



The T-dependence of Xi(qo;T) is only considered; that 
of X n (q a ;T) is ignored here. When the k dependence 
of </> 7 (k) is ignored, or £(k) given by Eq. (|2.48|) with 
7 (k) = 4> 1 is approximately used, it follows that 



^(Q,qo) 




-E— - 

N^k B T 

ha 



-sech 



2k B T 



(3.18) 



with 

esc(k) = V / e 2 (k) + (0 7 /0 7 ) 2 |iry r (k)e G (T)| 2 . (3.19) 

The function P(0,qo) is simply the Yosida function. 
When T = K, |e G (T)| > so that P(0,q ) = and 



X 1 (q o ;0K) = l/x s (0;0K). 



(3.20) 



Here, the T dependence of Xs(0) is explicitly shown. 
When T > T c but T < T K , on the other hand, e G (T) = 
so that 

P ^ * (t) - (3-21) and 



lea 



When the hybridization energy A(e) is constant, the local 
spin susceptibility at T — K is given by 



Xs(0) = 20 s p(O), 



(3.22) 



as is examined in Appendix [X] When T < Tk but no SC 
or no large pseudo-gap opens, i.e., T > 2T C , Eq. (|3.22p 
approximately holds, so that 

Xi(qo-,T)~(l-$ a /h)/x.(0-,T). (3.23) 

When Eq. (|2.55|) is satisfied, it follows that 

X 1 (q ;2T c )<l/x s (0;2T c ). (3.24) 

If the reduction of l/x s (0; OK) is large such that 

1/X S (0; OK) = Xtiqo; OK) < 1/ Xs (0; 2T C ), (3.25) 



the Knight shift is small; the T dependence of P(0, qo) 
is cancelled by that of l/x«(0;T). If the T dependence 
of l/x s (0;T) is small below T c or 2T C , the Knight shift 
is large, as it is in conventional superconductors. 

The Knight shift following the Yosida function is cer- 
tainly evidence that the symmetry of superconductivity 
is singlet. Since it can deviate from the Yosida function, 
however, a small Knight shift is never evidence that the 
symmetry of superconductivity is triplet if it occurs in 
the vicinity of the Mott M-I transition. 



IV. PHENOMENOLOGICAL THEORY 

In the vicinity of the Mott transition, two types of an- 
tiferromagnetic orderings are possible: screw or helical 
spin structure and sinusoidal spin structure or spin den- 
sity wave (SDW)i 72 ' 73 i 74 When the nesting of the Fermi 
surface of quasi-particles is significant in a metal, SDW 
is stabilized, which is presumably the case if antiferro- 
magnetism occurs in metallic cuprates. In general, SDW 
with wave number Q couples with charge density wave 
(CDW) with 2Q and lattice distortion with 2Q. Since the 
strong electron-phonon interaction arises from the modu- 
lation of the superexchange interactioni 41 ' 42 the coupling 
between SDW and lattice distortion must be strong but 
that between SDW and CDW is weak or is at least not 
strong. Then, Kondo-lattice theory should be extended 
to treat the coexistence of superconductivity, SDW, and 
lattice distortion. The extension is straightforward, as is 
discussed in Sec. [VlJ For the sake of simplicity, however, 
a phenomcnological theory is developed in this section. 

When an order parameter of SDW, that of lattice dis- 
tortion, and a coupling constant between them are de- 
noted by 



tuq = m_q = |mq|e 



U 2Q 



52Q 



-2Q = l U 2Q| 



'-2Q 



|& 2Q |e 4 



(4.1) 
(4.2) 



(4.3) 



with ±Q and ±2Q being their ordering wave numbers, 
a coupling term in Landau's free energy is given by 

= 2|& 2 QTOq u 2q| cos(6» 6 + 26 m 
The coupling term is minimal when 

cos(0 6 + 29 m - 6 U ) = -l. 



0. (4-4) 



(4.5) 



When an order parameter of superconductivity is de- 
noted by A S) Landau's free energy to be considered is 

F{T) = F + a s (T)\A s \ 2 + a rn (T) \m Q \ 2 



+C12Q|W2Q| -2|&2QmqU2Q| 



(4.6) 
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with 



and 



a s (T) = a s k B [T - T c ] , 



,(T) = a m k B [T - T m (Q)] 



(4.7) 



(4.8) 



Here, all the order parameters are defined to be di- 
mensionless. It is assumed that a2Q > 0, c s > 0, 
c m > |&2Q| 2 /a2Q, c sm > 0, a nd T c > T m (Q). According 
to the analysis in Sec. IIII A[ because of the reduction of 
T K by the opening of an SC gap, the coupling constant 
c sm can be much smaller in a superconductor in the vicin- 
ity of the Mott M-I transition than it is in a conventional 
one. Since the free energy (|4. 6|) is minimal when 



l"2Q 



°2Q | ,2 

\ m Q\ ■ 

0.2Q 



(4.9) 



it follows that 



F(T) = F +a s (T)\A s \ 2 + a m (T)\m Q \ 2 



^ S |A,| 4 



\ C Q \m Q \ 4 + c s 



with 

cq = c m - \b 2 Q\ 2 /a 2 Q > 0. 
Below T c , |A S | is non-zero. Define T^J(Q) by 



(4.10) 



(4.11) 



T£(Q) = T m (Q) 



[T c - T m (Q)] . 



(4.12) 

When T OT (Q) is large enough and c sm is small enough, it 
follows that TjJ(Q) > 0. In such a case, all of |A S |, |mq|, 
and |m2q| are non-zero below T^(Q). 
When T£(Q) < T < T c , it follows that 



A s \ 2 = -^k B [T c -T], 
2c s 



F(T) = F + Ft(T), 



F 1 (T) = -^k% [T c -T\ 
2c s 



The specific heat jump at T c is given by 

F,'2 



so that 



with 



ACi 



£1,2 rp 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



When < T < 2^(Q), it follows that 



cq 



k B [T* (Q) - T] , (4.17) 



^fc B [T c -T]-^|m Q | 2 , 



u 2q| is given by Eq. (|4.9p . and 

F(T)=F + J P 1 (r)+^ 2 (T;Q), 

with 



(4.18) 



(4.19) 



F 2 (T;Q) 



[am ^( Csm / C ^ . k %[T m (Q) -T} 2 . 



2 I C Q ~ C l m l C s] 



The specific heat jump at T^(Q) is given by 



AC 2 



[a m a s {c sm / c s )\ ^ 2 rp* 

2 / 
°Q c sm /c s 



(4.20) 



(4.21) 



When T m (Q) is mainly determined by J s (q) defined 
by Eq. (|2.4p . T m (Q) is maximal at wave numbers in the 
vicinity of 



Qm = (±7r/a,±7r/o) 



(4.22) 



When the enhancement of T m (Q) by the strong electron- 
phonon interaction is large, T m (Q) can be maximal at 
wave numbers different from Qa/- Denote one of the 
wave numbers by Qo- The Neel temperature is given by 



Tn — ?n(Qc 



(4.23) 



The ordering wave number Q below Tn is determined by 
minimizing F(T) or F 2 (T; Q) as a function of Q. When 
only the Q dependences of T m (Q), a2Q, and |&2q| are 
considered, it is enough to minimize 



P 1 * (Q) - Tf 



Cm - |&2Q| 2 / a 2Q - C 2 sm /C s 



(4.24) 



instead of F 2 (T;Q). When the electron-phonon inter- 
action is strong, therefore, the ordering wave number of 
SDW can be different from Qm at Tn and it can change 
with decreasing T below T^. 



V. CHECKERBOARD AND 
ZERO-TEMPERATURE PSEUDO-GAP IN 
CUPRATE SUPERCONDUCTORS 

The DOS depends on the position r in the checker- 
board phase of cuprate superconductors, whose period is 
of about 4a x 4a 49 i 50 1 5i 1 52 1 53 i 54 T h e DOS defined by 



p(£;Qc)= (2^ / dreJQc >( £ ; r ) 



(5.1) 



with Q c ~ (±7r/2a,0) or (0, ±7r/2a) being one of wave 
numbers of the checkerboard and p(e; r) being DOS at 
r summed over spins, is almost symmetric with respect 
to e = such that p(e; Q c ) — p(— e;Q c )- 51 It cannot 
be explained in terms of CDW by itself or CDW in- 
duced by lattice distortion, for which p(e; Q c ) is almost 
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asymmetric with respect to e — such that p(e; Q c ) ~ 
— p{— e; Q c ). The ZTPG seems to be always accompanied 
by the checkerboard , 53 ! 55 ^ which implies that the ZTPG 
and the checkerboard below T c are different aspects of 
a phase. Then, a possible interpretation for them is a 
double- Q CDW or SDW in a conventional or exotic su- 
perconducting state: Cooper-pair density wave, a Wigner 
solid of Cooper pairs, and so n»fflfflffl2ffli 
As is studied in Sees. IIII Al and HVl SDW can coexist with 
^7-wave superconductivity. It is shown that the almost 
symmetric p(e; Q c ) can be explained if CDW with an in- 
commensurate Q c is induced by SDW with iQ c + G, 
G being a reciprocal lattice vector, as a second-order ef- 
fect of the SDWj^ It is also shown that a similar fine 
structure to that of ZTPG can be reproduced if the co- 
existence of superconductivity and SDW is assumed. 84 
The purpose of this section is to propose that, provided 
that T <T C . the checkerboard and the ZTPG phases are 
different aspects of a coexistence phase of cfry-wave super- 
conductivity, a double-Q SDW, a double-Q CDW, and a 
double- Q lattice distortion. 

Denote ordering wave numbers of SDW by 

Q = Qm + S{tt/ci), (5.2a) 

with Qm = (±7r/a, ±7r/a) and 

S = (±S X , ±S y ), (5.2b) 

and those of the CDW and the lattice distortion by 

2Q = <*(27r/a). (5.2c) 

According to previous papers ; 41 ! 42 it is straightforward to 
show that |&2q| has a peculiar dependence on 2Q or 5. 
The mode with 2Q = Qjy./ is the breathing mode. Since 
it does not modulate the superexchange interaction, |&2q| 
are small around 2Q = Q M or S = (±1/2, ±1/2). The 
mode with 2Q = Qx, with Qx = (±7r/a, 0) or (0, ±n/a) 
is the half-breathing mode. Since it modulates strongly 
the superexchange interaction, &2Q is maximal around 
2Q = Qx, or S = (±1/2,0) and S = (0,±l/2)^ Then, 
it is probable that the enhanced T m (Q) by the electron- 
phonon interaction is maximal on lines of 

< |<J X | < 1/2, 5 y = 0, (5.3a) 

and 

S x = 0, < 1^1 < 1/2, (5.3b) 

in the plane of (S x , 5 y ) and the free energy (|4.24p is mini- 
mal on the lines. It is possible that the free energy (|4.24|) 
is minimal around midpoints of the lines such as 

*i = (±1/4,0), (5.4a) 

and 

6 2 = (0, ±1/4) . (5.4b) 



Two SDWs with Qi = Q M + S^n/a) and Q 2 = 
Qm + ^{jr/a) are equivalent to each other in the square 
lattice. According to Kondo-lattice theory of SDW^i a 
double-Q SDW with Qi and Q2 must be stabilized rather 
than a single-Q SDW with Qi or Q 2 since an additional 
stabilization energy appears when magnetization of dif- 
ferent Q components are orthogonal to each other in the 
double-Q SDW. The double-Q SDW with Qi and Q 2 is 
inevitably accompanied by double-2Q lattice distortion 
and CDW with 2Qi and 2Q 2 . The period of the double- 
2Q lattice distortion and CDW, which is of about 4ax4a, 
is the observed period of the checkerboard structure 

According to the previous paper^s 2 . phonons with 
2Qi = di(2ir/a) and 2Q 2 = <52(27r/a), which are of the 
Cu-0 bond-stretching mode, can become soft provided 
that antiferromagnetic spin fluctuations with Qi and 
Q2 are developed. The softening of the bond-stretching 
mode is actually observedi 36 ! 37 In reverse, the softening 
enhances the instability of the double- Q SDW. 

If order parameters of d-y-w&ve superconductivity and 
the double-Q SDW coexist homogeneously with each 
other, the dispersion relation of quasi-particles is cooper- 
atively renormalized by the two order parameters or the 
opening of two types of gaps. Then, DOS can have a fine 
gap structured The fine structure must be ZTPG. 

VI. DISCUSSION 

It is straightforward to extend Kondo-lattice theory to 
study the coexistence of superconductivity, SDW, lattice 
distortion, and CDW. When the wave number Q of SDW 
is incommensurate, the Green function is given by 

f 1/ksT ■ I - «+ \ 

G a (ie n , k) = - jf dre^ T (T T A k(r {0), 4») , (6.1) 

with A k(T being the Hermitian conjugate of 

^kcr = (' ' ' ' a k-2Qcr' fl k-Qcr' a k<r' a k:+Qer' fl k+2Qcr' 

1 a-k-2Q-er, a-k-Q-er, O-k-er, 

Q-k+Q-tT, a-k+2Q- CT , ' ' ' )• (6.2) 

The self-energy and order parameters appear in the di- 
agonal and off-diagonal parts, respectively, of Eq. (|6.1[) . 
The formulation is almost in parallel with that in Sec.HIl 
Since the Fock-type term or the RVB mechanism, 
which suppresses the appearance of SDW, is reduced 
when the life-time width of quasi-particles is large, SDW 
more easily appears in more disordered systems than it 
does in less disordered systems^ It is possible that dis- 
order is much larger on surface Cu02 layers than it is 
on deep Cu02 layers. In such a case, the checkerboard 
and ZTPG phase may exist only on disordered surface 
Cu0 2 layers and it appears inhomogencously on them. 
Experimentally, it is certain that the checkerboard and 
ZTPG phase appears on surface Cu02 layers. It should 
be examined whether it can appear on deep Cu02 layers. 
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The band center of LHB is around e ~ — W/2. Ac- 
cording to Hubbard's theory^ p(e ~ —W/2) ~ l/W, 
so that p(e ~ —W/2) ~ 0.3 eV _1 for cuprates. When 
ImS CT (E + iO) and Rei? (T (e + iO) are ignored in Eq. (|2.17[) . 
A(e ~ — W/2) ~ 1 eV. When they are considered, 
it is plausible that A(e ~ —W/2) ~ | x 1 eV, i.e., 
Ai ~ 300 meV. On the other hand, A(0) given by 
Eq. (|2.19[) decreases as p(0) decreases provided that 
Rei? (T (+i0)| > np(0). When a pseudo-gap structure de- 
velops in p(e) above T c , A(e) must also have a pseudo-gap 
structure. Experimentally, a half of the pseudo-gap is as 
large as eo — 50 meV. Then, it is plausible that 

eo/Ai ~ 0.2, (6.3) 

in the pseudo-gap phase. According to the results shown 
in Fig. [TJ it is possible that a double-Q SDW or the 
checkerboard appear at T ~ T c or T > T c if the pseudo- 
gap is so largely developed that not only Eq. (|6.3[) but 
also Ao/Ai <C 1 are satisfied. The checkerboard is actu- 
ally observed in the pseudo-gap phase above T c 4& 

Since two modes of antiferromagnetic spin-fluctuations 
with qi and q2 inevitably couple with a phonon mode 
with qi + q2 due to the electron-phonon interaction, spin 
excitation spectra are renormalized by phonons. It is in- 
teresting to examine the relevance of such a spin-phonon 
coupled mode to the so called resonance mode, which is 
observed in neutron inelastic scattering experiment 42 It 
is also interesting to examine whether the Knight shift 
follows the Yosida function in under-doped cuprate su- 
perconductors, in which the pseudo-gap, the checker- 
board, or ZTPG is observed. 

When the Brillouin zone (BZ) is folded by SDW, CDW, 
and lattice distortion, Cooper pairs around one of the BZ 
edges or Cooper pairs with non-zero total momentum can 
coexist with Cooper pairs around the BZ center 4^ A co- 
existence phase of such two types of Cooper pairs is an- 
other type of Cooper-pair density wave. It is interesting 
to search for such a coexistence phase. 

VII. CONCLUSION 

According to Kondo-lattice theory of high-temperature 
superconductivity in cuprate superconductors, which is 
developed in a series of previous papers, the normal state 
is the Fermi liquid, which is mainly stabilized by the 
Kondo effect and the resonating valence bond (RVB) 
mechanism, and cfy-wave superconductivity occurs on 
Cu02 planes mainly due to the superexchange interac- 
tion between-nearest neighbor Cu ions. In this paper, 
the Kondo-lattice theory is extended to treat properly 
the Kondo effect or the quenching effect of magnetic mo- 
ments by single-site quantum spin fluctuations in a d'y- 
wave superconducting state. It is shown that, since the 
quenching effect is weakened by the opening of a dj- wave 
superconducting gap, spin density wave (SDW) can ap- 
pear in the superconducting state and the Knight shift 
can deviate from the Yosida function to be small. 



On the basis of the Kondo-lattice theory, together with 
the previous theory on the electron-phonon interaction 
on Cu02 planes that arises from the modulation of the 
superexchange interaction by phonons, a phcnomcnolog- 
ical theory is developed to study the coexistence of su- 
perconductivity, SDW, and lattice distortion. Because 
of the strong electron-phonon interaction, the lattice dis- 
tortion inevitably appears and the wave numbers of the 
SDW can be different from (±7r/a, ±7r/a), with a being 
the lattice constant of Cu02 planes, for which the su- 
perexchange interaction is maximal. 

On the basis of the phenomenological theory, it is pro- 
posed that the checkerboard and the zero-temperature 
pseudo-gap (ZTPG) in the superconducting phase must 
be different aspects of a coexistence phase of G?7-wave su- 
perconductivity, a double-Q SDW, a double-2Q lattice 
distortion, and a double-2Q charge density wave (CDW). 
The checkerboard can be explained if the two wave 
numbers of the SDW are different from (±7r/a, ±7r/o) 
because of the strong electron-phonon interaction, i.e., 
those of the SDW are Qi ~ (±3n/4a, ±n/a) and Q 2 ~ 
(±7r/o, ±37r/4a) and those of the lattice distortion and 
the CDW are 2Qi ~ (±7r/2a,0) and 2Q 2 ~ (0,±7r/2a). 
The fine structure of the ZTPG must be mainly due 
to the opening of superconducting and and SDW gaps. 
On the basis of previous theories, it is discussed that 
the stabilization mechanism of the double-Q structure 
or the checkerboard rather than a single-Q structure or 
the stripe is that an additional stabilization energy arises 
in the double-Q SDW when magnetic moments of differ- 
ent Q waves are orthogonal to each other, the double-2Q 
CDW must be mainly induced by the double-Q SDW 
rather than the double-2Q lattice distortion since the 
strong electron-phonon interaction, which arises from 
spin channels, can play no significant role in the induc- 
tion of the CDW, and the strong electron-phonon inter- 
action can play no significant role either in the occurrence 
of d'y-w&ve superconductivity itself since the phonon- 
mediated attractive interaction is strong between next 
nearest neighbors but is weak between nearest neighbors. 
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APPENDIX A: SUSCEPTIBILITY OF THE 
ANDERSON MODEL 

Consider the Anderson model (|2. 12[) at T = K in the 
presence of an infinitesimally small Zeeman energy —ah. 
The Green function for d electrons is given by 



ie n + p, - e d + <rh- t, a (ie n ) - L a (ie n ) ' 

(Al) 



12 



and that for conduction electrons is given by 

G CCT (ie„;k,k') = g ca {ie n , k)5 kk / + g ca {ie n , k) 

xVkG^eJV^^k'), (A2) 



with 



and 



# C(T (i£ n ,k) 



ie n + jl- E c (k) ' 



^<r(i£n) = i |Vk| 2 g CCT («£„,k). 



The number of d electrons is given by 



del lmG a (e + i0), 



(A3) 



(A4) 



(A5) 



and the change in the number of conduction electrons 
due to the hybridization is given by 



An r „ = 



de—lm 

7T 



ae 



(A6) 



The spin susceptibility due to d electrons is defined by 
Xs{ti) = (d/dh)Y J ™d„, (A7) 

and that due to conduction electrons is defined by 

Axs(0) = (d/dh)J2^ri ca . (A8) 



It is straightforward to prove that 



X.(0) + Ax.(0) = 2(M0), 



(A9) 



with cj) s the expansion coefficient defined by Eq. (|2.20|) 
and p(0) given by Eq. (|2~Td]1 with Eq. f2TT5]) , 

When the SC gap opens in the Hubbard or t- J model, 
p(0) = in the mapped Anderson model so that 



X s (0) + Ax s (0) = 0. 



(A10) 



On the other hand, X s (0) is enhanced, as is studied in 
Sec. IIII Al Then, A\ s (0) is also enhanced in such a way 
that polarizations of conduction electrons exactly cancel 
those of d electrons. Anderson's compensation theorem 
does not hold, i.e., Ax s (0) ^ in this case. 

Consider a model where the hybridization energy de- 
fined by Eq. (j2 . 13[) is given by 



A(e) = A(0)r 2 /(e 2 + r 2 ), 
or L a {e + iO) = A(0)r/(e + iT). It follows that 



(AH) 



8L a {e + iQ) 
de 



-A(0) 



{e + iYf 



(A12) 



In the limit of V — > +oo, A(e) is constant. In such a case, 
Eq. (|A"l2|) vanishes so that An ca = 0. Then, Ax s (0) = 
or Anderson's compensation theorem holds. 
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